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Abstract 

We investigate a possible scheme for entangling two mode thermal fields 
through the quantum erasing process, in which an atom is coupled with 
two mode fields via the interaction governed by the two-mode two-photon 
Jaynes-Cummings model. The influence of phase decoherence on the entan- 
glement of two mode fields is discussed. It is found that quantum erasing 
process can transfer part of entanglement between the atom and fields to 
two mode fields initially in the thermal states. The entanglement achieved 
by fields heavily depends on their initial temperature and the detuning. 
The entanglement of stationary state is also investigated. 
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I. INTRODUCTION 



Entanglement, an important resource for quantum information processing [1], is one 
of the most prominent nonclassical properties in quantum theory. Entanglement can 
exhibit a nonlocal correlation between quantum systems that have no classical inter- 
pretation. Recently, much attention has been focused on the entanglement in bipartite 
or multipartite systems in which the subsystems are initially in thermal equilibrium 
[2,3,4,5]. Arnesen et al. have shown that a natural entanglement arises in Heisenberg 
spin chain in thermal equilibrium, and the entanglement can be improved by increasing 
the temperature [2]. Instead of attempting to shield the system from the environmental 
noise, Plenio and Huelge [3] use white noise to play a constructive role and generate 
the controllable entanglement by incoherent sources. Similar work on this aspect has 
also been considered by other authors [4,5]. However, very little attention has been 
paid to the study of entangling two mode thermal fields. In this paper, we investigate 
a possible scheme for entangling two mode thermal fields through the quantum erasing 
process, in which an atom is coupled with two mode fields via the interaction governed 
by the two-mode two-photon Jaynes-Cummings model. The influence of phase deco- 
herence on the entanglement of two mode fields is discussed. It is found that quantum 
erasing process can transfer part of entanglement between the atom and fields to two 
mode fields initially in the thermal states. The entanglement achieved by fields heavily 
depends on their initial temperature and the detuning. The term "quantum eraser" 
[6] was invented to describe the loss or gain of interference or, more generally quantum 
information, in a subensemble, based on the measurement outcomes of two comple- 
mentary observables. It was reported that the implementation of two- and three-spin 
quantum eraser using nuclear magnetic resonance, and shown that quantum erasers 
provide a means of manipulating quantum entanglement [7]. The quantum erasing 
process discussed in this paper is implemented by measuring the polarizing vector of a 
two-level atom coupling with two mode quantum fields. The project measurement of 
an atom has been extensively studied both in the theoretical and experimental aspects. 
This paper is organized as follows. In Sec. II, we study the system in which an atom is 
coupled with two mode fields via the interaction governed by the two-mode two-photon 
Jaynes-Cummings model by making use of the dynamical algebraical method [8,9] and 
find the exact solution of the master equation for the system with phase decoherence. 
Based on the exact solution, we then propose a possible way to entangle two mode 
thermal fields through the quantum erasing process, which is realized by measuring 
the atom. In Sec. Ill, we use the log-negativity to characterize the entanglement be- 
tween two mode fields. It is shown that quantum erasing process can transfer part of 
entanglement between the atom and fields to two mode fields initially in the thermal 
states. The entanglement achieved by fields heavily depends on their initial tempera- 
ture and the detuning. A conclusion is given in Sec. IV. 
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II. SOLUTION OF AN ATOM COUPLES TO TWO 



THERMAL FIELDS WITH PHASE DECOHERENCE 

We consider the two-mode two-photon Jaynes-Cummings model [10]. The Hamiltonian 
for the model can be described by (h — 1), 

H = c^iajai + uj 2 a 2 a 2 + ^a z + g(a 1 a 2 a + + a\a\<jJ), (1) 

where a z and o± are the atomic spin flip operators characterizing the effective two-level 
atom with transition frequency u and a± (a 2 ), a\ (a 2 ) are annihilation and creation op- 
erators of the first (second) mode light field of frequencies ui (oo 2 ) respectively. The 
Hamiltonian (1) ignores Stark shifts and the parameter g is the atom-field coupling 
constant. 

It is easy to see that there exist two constants of motion in the Hamiltonian (1), 

K l = a[ax H — , K 2 = a 2 a 2 H — , (2) 

which commute not only with Hamiltonian but also with operators aia 2 a + and a\a 2 cr-. 
We can introduce the following operators 

S =— S = aia2<T+ Q = a i a 2°"- (3) 
° 2 ' + ^/KJT 2 ~ 

The operators S± and S satisfy the following commutation relations 

[S ,S ± ] = ±S ± , [S + ,S_]=2S , (4) 

where S± and S are the generators of the su(2) algebra. In terms of the su(2) gener- 
ators, we can rewrite the Hamiltonian (1) as 

H = W\(Ki - 1) + uj 2 (K 2 - l -) + AS + g^K 1 K 2 (S + + SI), (5) 

where A = uj — oj\ — uj 2 . With the help of the su(2) dynamical algebraic structure, we 
can diagonalize the Hamiltonian (5) by introducing a unitary transformation 

C/ = exp[^%^(S + -S_)| (6) 



with 8(K 1 ,K 2 ) = aictan(2gy/K 1 K 2 / A) , and get transformed Hamiltonian 

H' = UHU^ = oj^K, - 1) + lo 2 (K 2 - 1) + 2fi(^i, K 2 )S , (7) 



where tt(K u K 2 ) = ^A 2 /4 + g 2 K x K 2 . 

In this paper, we consider the phase decoherence mechanism only. In this situation, 
the master equation governing the time evolution for the system under the Markovian 
approximation is given by [11] 

% = -i[H,p]-l[H,[H,p]], (8) 



where 7 is the phase decoherence coefficient. Noted that the equation with the similar 
form has been proposed to describing the intrinsic decoherence [12]. The formal solution 
of the master equation (8) can be expressed as follows [13], 

P(t) = T, {S TTM k (t)p(0)M^(t), (9) 

k=0 

where p(0) is the density operators of the initial atom-field system and M k (t) is defined 
by 

M k (t) = H k exp{-iHt) exp(-^# 2 ). (10) 

By means of the SU(2) dynamical algebraic structure, we obtain the explicit expression 
for the operator M k 

M k (t) = U^H' k exp(-iH't)exp(-^H >2 )U 
= \[fl exp(-i/ + f) exp(-^t) + f k _ eM-if-t) exp(-^)] 

(ii) 

where f± = uj\{Ki — 1/2) + uj 2 {K 2 — 1/2) ± Q(Ki,K 2 ). Firstly, we assume that the 
cavity fields are initially in two-mode Fock states \n1n2), and the atom is in the excited 
state |e). The time evolution of p(t) can be written as follows, 

= it 2 + + ( 2 " J^) e ~ 2imll ' n2 cos2Q niin2 t]\n 1: n 2 )(n 1: n 2 \ <g> \e)(e\ 

4 Zi 'ni,n 2 Zi 'ni,n 2 

1 q 2 (ni + l)(n 2 + 1) r % ,, n 2 „ ,, . , , . , , 

+ j 02 L[2-2e- 2 ^^ cos2fi ni , n2 t]|n 1 + l,n 2 + l)(n 1 + 1, n 2 +l\®\g){g\ 



^yTn^flK^Tl) A _ 27tn2 oo 



4^ni,n2 ^ni,c- 



+ 2ie 27 *^i.«2 sin2n f}|ni,n 2 )(ni + l,n 2 + 1| <g> \e)(g\ 



+ ^"': 1K "' + " l ^[l - e-<,™ co.2lW] 
- 2ie- 27tn "i." 2 sin2fi„ lin2 t}|n 1 + l,n 2 + l)(m,n 2 | ® \g){e\, (12) 



where fi nii „ 2 = ^/A 2 /4 + g 2 ( ni + l)(n 2 + 1). In the basis {|1, 1) = |m, n 2 )(g)|e), |0, 1) = 
|ni + l,n 2 + 1) <g> |e),|l,0) = |n 1? n 2 ) <g> |#),|0,0) = \m + l,n 2 + 1) <g> |#),}, p(t) 
can be regarded as a two qubit mixed state. Then, a quantum erasing is applied to 
this system by making a project measurement of the atom on the basis (cos||e) + 
e 1 ^ sin fig), cos f |g) — e~ % ^ sin f |e)}. It is easy to verify that two fields will get the same 
amount of entanglement corresponding to two different measurement outcomes if the 
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value of 9 is ir/2. In this case, both the probabilities of two projection measurement 
results are \. So we can only consider the entanglement of one of the projection results 
instead of average entanglement between the two fields after the measurement. If the 
measurement projects the state of the atom onto cos||e) + e"^sin||g), the residual 
state of two mode fields is expressed by (unnormalized) 

1 9 A 2 A 2 2 
Pf (n l7 n 2 , t) = - cos 2 - [2 + — + (2- — ) e - 2 ^U 2 cos 2Q nun2 t] | \n u n 2 ) (m, n 2 \ 

1 9 9 q 2 (ni + l)(no + 1) r 9/v . 2 ^ 

+ - sin 2 - y V A ^[2 - 2e- a ^»i.»2 cos 2fi nii „ 2 t] | |m + 1, n 2 + 1) (m + 1, n 2 + 1| 

4 ^ ni,«2 



+ 2ie" 27 * Q "i."2 sin2^ 1)n2 t}||ni,n 2 )(ni + l,n 2 + 1| 



- 2ie- 27 * n »i.»3 sin2ft niina *}|K + l,n 2 + l)(n 1? n 2 | (13) 
For the initial two mode thermal fields, the output state of two fields is replaced by 



mi l m 2 2 



"* (() = tf Ju (l + r„ 1 )«.«(l + ^)-« P/(n " ' 12 ' 



(14) 



where Af = {E^, n2 =o (i+mi)™+ 1 aW" 2+l Tr M Wl ' n 2, ^j}' 1 is the normalization con- 
stant, and rhi = [exp(fau)i) — (i = 1, 2) is the mean photon number of the ith mode 
thermal field at the inverse temperature fa. 



III. THE LOG-NEGATIVITY OF TWO MODE FIELDS 

In order to quantify the degree of entanglement, we adopt the log-negativity N(p) to 
calculate the entanglement between two mode fields, which is defined as [14] 

N(p) =log 2 ||p r ||, (15) 

where p r is the partial transpose of p and ||p r || denotes the trace norm of p r , which is 
the sum of the singular values of p T . 



For the unnormalized density operator in Eq.(13), it is easy to derive its stationary 

log- negativity which is given by log 2 [l + 2| sin 9 9 ^ <yni + a ^Tote 1 1 • ^ or simplicity, we will 
set the value of 9 as | throughout the following calculation. First of all, one important 
fact should be pointed that the entanglement between two initial thermal fields can not 
arise if the quantum erasing processing is not applied and the degree of freedom of the 
atom is simply traced. The entanglement between the fields is partly transferred from 
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Figure 1: The stationary state log-negativity N(pf) of the density operator p/(oo) is 
plotted as a function of the mean photon number fhi = fh 2 = a of initial thermal fields 
and the detuning A with g = 0.5 and 9 = | . 

the entanglement between the atom and the fields through the quantum erasing. So, 
all of the following discussions concerning the entanglement between two fields at any 
time t are based on the presumption that a projection measurement is just acted on the 
atom at the time t~. In Fig.l, the stationary state log-negativity N(pf) of the density 
operator p/(oo) is plotted as a function of the mean photon number fh\ = m 2 = ol 
of initial thermal fields and the detuning A. Fig.l shows that, in the resonant case 
A = 0, there is not any entanglement in the stationary state. In the off-resonant case, 
the entanglement decreases with a, and eventually disappears as the value of a goes 
beyond a threshold value which is dependent on the detuning. A natural question 
will arise how the entanglement behaves when one mode is initially in the vacuum 
state and the other mode is in thermal state. In Ref. [15], the authors indicated that 
the subsystem purity can enforce the entanglement. So, it is easy to understand the 
result displayed in Fig. 2, where the stationary state log-negativity N(pf) of the density 
operator p/(oo) is plotted as a function of the mean photon number fh 2 = a of initial 
thermal field of the second mode and the detuning A with fhi = 0, i.e., the first mode 
is initially in a vacuum state. We can find that the entanglement always exists for 
any high temperature of the second mode in the off-resonant situation. One may also 
conjecture that the stationary state entanglement can increase with the difference of 
the mean photon numbers of two thermal fields as the value of fhi + fh 2 is fixed. This 
seems to be true, and can be seen from the Fig. 3, in which we depict the stationary 
entanglement as a function of fh\ and fh 2 with g = 0.5 and A = 1. The stationary 
entanglement always increases with the value of \fhi —fh 2 \ along any line characterized 
by fhi + fh 2 = const.. We conjecture this phenomenon exists in a wide class of systems, 
including the thermal modes in different thermal reserviors effectively coupled by the 
qubits in a quantum register. 

In the resonant case, there is not any stationary state entanglement between the 
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Figure 2: The stationary state log-negativity N(pf) of the density operator p/(oo) is 
plotted as a function of the mean photon number m 2 = a of second mode thermal field 
and the detuning A with g = 0.5, fh\ = and 9 = f . 




Figure 3: The stationary state log-negativity N(pf) of the density operator p/(oo) is 
plotted as a function of the mean photon number rhi and rfi2 of initial thermal fields 
with g = 0.5, A = 1 and 6 = f . 
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Figure 4: The log-negativity N(pf) of the time evolution density operator Pf(t) is 
plotted as a function of the mean photon number fhi = m 2 = ot of initial thermal fields 
and the time t with g = 0.5, 9 — |, 7 = 0.5 and A = 0. 

two modes. Nevertheless, the entanglement still arise in the forepart of the evolution, if 
either the initial temperature of the thermal fields or the phase decoherence coefficient 
are not too large. In Fig.4, the log-negativity N(pj) of the time evolution density 
operator Pf{t) is plotted as a function of the mean photon number rhi = m 2 — ot 
of initial thermal fields and the time t. It is shown that the two-mode fields can get 
entangled in the beginning of the time evolution, and become disentangled due to the 
presence of decoherence. However, in the off-resonant case, the entanglement is robust 
against the phase decoherence. Fig. 5 clearly displays how the two initial thermal 
fields get entangled and eventually evolve into a stationary entangled state. When 
the two fields are initially in thermal states, the higher the temperature, the later 
the onset of entanglement between two fields. There will not be any entanglement 
appearing between two fields as their initial temperature exceeds certain threshold 
value which depends on the decoherence coefficient, the coupling strength and the 
detuning. Furthermore, we plot the log-negativity as the function of the time and the 
mean number difference 5 = \fhi— fri2\ with a fixed mean number sum rhi+rhz = 1. It is 
shown that the log-negativity increases with the number difference at any time. In the 
case with uj\ = 002, increasing the difference of initial temperature of two thermal field 
results in enlarging the value of 5. So, one can improve the entanglement by increasing 
the temperature difference in the situation that the total energy of the initial thermal 
fields is fixed. The novel phenomena that increasing the temperature difference of 
the thermal fields will improve their entanglement may have some applications in the 
quantum information processing, in which some subsystems are initially in thermal 
equilibrium. 



8 





Figure 5: The log-negativity N(pj) of the time evolution density operator is 
plotted as a function of the mean photon number rhi = m 2 = a of initial thermal fields 
and the time t with g = 0.5, 9 — |, 7 = 0.5 and A = 1. 




Figure 6: The log-negativity N(pf) of the time evolution density operator is 
plotted as a function of the mean photon number difference 5 = \rhi — 777.2 1 of initial 
thermal fields and the time t with g = 0.5, 9 — ^, 7 = 0.5, 777! + 77i 2 = 1 and A = 1. 
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IV. CONCLUSION 



In this paper, we investigate a possible scheme for entangling two mode ther- 
mal fields through the quantum erasing process, in which an atom is coupled with 
two mode fields via the interaction governed by the two-mode two-photon Jaynes- 
Cummings model. The influence of phase decoherence on the entanglement of two 
mode fields is discussed. It is found that quantum erasing process can transfer part of 
entanglement between the atom and fields to two mode fields initially in the thermal 
states. The entanglement achieved by fields heavily depends on their initial tempera- 
ture and the detuning. The entanglement of stationary state is also investigated. It 
is interesting to study the entanglement in a similar scheme in which the two-mode 
two-photon Jaynes-Cummings model is replaced by the two-mode Raman coupling 
Jaynes-Cummings model. Both schemes can be easily realized in the two-dimensional 
ion trap. The details will be discussed elsewhere. 
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